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We study the self force acting on a particle endowed with scalar charge, which is held static 
(with respect to an undragged, static observer at infinity) outside a stationary, axially-symmetric 
black hole. We find that the acceleration due to the self force is in the same direction as the 
black hole's spin, and diverges when the particle approaches the outer boundary of the black hole's 
ergosphere. This acceleration diverges more rapidly approaching the ergosphere's boundary than 
the particle's acceleration in the absence of the self force. At the leading order this self force is a 
(post)^-Newtonian effect. For scalar charges with high charge-to-mass ratio, the acceleration due 
to the self force starts dominating over the regular acceleration already far from the black hole. 
The self force is proportional to the rate at which the black hole's rotational energy is dissipated. 
This self force is local (i.e., only the Abraham-Lorentz-Dirac force and the local coupling to Ricci 
' curvature contribute to it). The non-local, tail part of the self force is zero. 
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I. INTRODUCTION AND SUMMARY 



Recently there has been much interest in the calculation of self interaction of particles in curved spacetime. This 
growing interest is motivated by the prospects of detection of low-frequency gravitational waves in the not-so-distant 
^ \ future by space-borne gravitational wave detectors such as LISA. The main challenge is to compute the orbital 
evolution, and the resulting wave forms, from a compact object orbiting a supermassive central black hole in the 
extreme mass ratio case. The motion of the compact object, which may be construed as a structureless particle, is 
, geodesic in the limit of zero mass. However, when the particle is endowed with a finite (albeit small) mass, its motion 
' is changed. Spacetime is now determined by the energy-momenta of both the black hole and the particle, and the 
T-H , latter follows then a geodesic of the new spacetime, which is perturbed by its own energy-momentum relative to the 
original spacetime. An alternative viewpoint is to consider the motion of the particle as an accelerated, non-geodesic 
motion in the unperturbed spacetime of the central black hole. Whereas this latter approach is less in the spirit of 
Einstein's General Relativity, which "eliminated" gravitational forces in favor of geometry, than the former approach, 
it has the advantage that the unperturbed spacetime is often very simple. (For many interesting cases it is, e.g., the 
(^JT), stationary and axially-symmetric spacetime of a Kerr black hole.) We are thus led naturally to translate the problem 
^ ' of finding the orbital evolution of the particle to the following question: what are all the momentary forces which act 
• on the compact object? (In the absence of external forces, there would be just the self force.) Obviously, knowledge 

■ of all the forces which act on an object in a given spacetime allows for the computation of the orbit and consequently 
J-j \ also for the computation of the emitted gravitational waves. 
. - - 1 The problem of finding the self force which acts on a particle in curved spacetime is not easy. The reason is that 
in curved spacetime, due to the failure of the Huygens principle, the retarded Green's function associated with the 
particle's field has support inside the future light cone, and in particular also on the future world line. (The physical 
origin for this phenomenon is the scattering of the emitted waves off the spacetime curvature.) The part of the Green's 
function inside the light cone has been dubbed its "tail" part, and it is the calculation of this tail part of the retarded 
Green's function which is the greatest challenge in the computation of the self force (a.k.a. the radiation reaction 
force). A difficulty associated with the calculation of the tail part of the Green's function is the prescription used 
in order to separate the tail part, which affects the motion of the particle, from the instantaneous, divergent part, 
which arises from the typical divergence of the particle's field in the coincidence limit of the source for the field and 
the evaluation point. 

Several approaches have been proposed for the regularization of the self force (SF). One approach was first suggested 
by Dirac for a pointlike electric charge in flat spacetime 0], and later used by DeWitt and Brehme for the case of an 
electric charge in curved spacetime [|| and by Mino, Sasaki, and Tanaka for the case of a pointlike particle coupled 
to linearized gravity The idea is to impose local energy and momentum conservation on a tube surrounding the 
particle's world line, and to integrate the conservation laws across the tube, thus obtaining the particle's equations 
of motion, including the SF effects. The divergent piece of the self force is then removed by a mass-regularization 
procedure. A second, axiomatic approach, which leads to the same expression as the approach described above, was 
proposed by Quinn and Wald for electromagnetic and gravitational self forces B and by Quinn for the scalar field case 
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1^. According to this approach, the regularization of the SF is performed by comparing the forces in two different 
spacetimes. 

For the case of a particle coupled to a (minimally-coupled, massless) scalar field, the total SF which acts on the 
particle is given by 

i^SF = 1 2 _ y^y^^o) + 1 27^^ + q2 r""' d^'v^G-*[z"(r), z"(r')] , (1) 

6 D £^0 J-00 

where 7?.^ = -R^qM" + Rapu'^u^u^j, — ^Ru^. Here, is the Ricci tensor, R is the curvature scalar, and u" is the 
particle's four-velocity. An overdot denotes (covariant) differentiation with respect to proper time r, q is the particle's 
scalar charge, G"''^*[z"(t), 2;"(r')] is the retarded Green's function, and z"(r) is the particle's world line. The total 
force which acts on the particle is the sum of external forces (e.g., forces which result from external scalar fields) and 
the SF (|l|). The SF in Eq. (|^) has three contributions: the first is a local, Abraham-Lorentz-Dirac (ALD) type force. 
The ALD force consists of two terms, a term proportional to the proper time derivative of the four-acceleration, called 
the "Schott part" of the ALD force, and a term proportional to the four-acceleration squared, which we shall call here 
the "damping part." The second term of Eq. (|l|) comes from local coupling of the particle to Ricci curvature. This 
term preserves the conformal invariance of the SF. The two local terms, namely the ALD and Ricci-coupled terms, 
constitute the local part of the SF, Fj^'^^^. The third term in the SF is the non-local tail term, which involves an 
integration over the gradient of the retarded Green's function along the entire past world line of the particle. 

A number of practical methods have been used in order to find the effects of the tail part of the SF. One such 
approach avoids the calculation of the SF by applying balance arguments to quantities which are constants of motion 
in the absence of radiation reaction (i.e., the energy and angular momentum). This method has been successful to 
quite high relativistic order when the black hole is spherically symmetric, or when the motion is circular or equatorial 
around a spinning black hole, because in these cases the orbital evolution is uniquely determined by the rate of 
change of energy and angular momentum, and global conservation laws imply that the fluxes of energy and angular 
momentum across a distant sphere and down the event horizon of the black hole equal the rate of loss of energy and 
angular momentum by the particle, respectively. However, this method is useless for cases where there is non-trivial 
evolution of Carter's constant, which is a non-additive constant of motion in the absence of radiation reaction. This 
method has in addition two other problems. First, it involves averaging over the orbital period. Thus, when the 
orbital evolution is fast and the time scale for the orbital evolution becomes comparable with the time scale for the 
orbital period, this method becomes inaccurate. Second, this method is sensitive to the slowest decaying part of the 
field at large distances. Consequently, it handles well the dissipative, radiative part of the field, but completely ignores 
the conservative part of the SF. Although the latter does not cause any net loss of energy and angular momentum, it 
is still important for the orbital evolution |^]. In fact, even in the absence of dissipation, the conservative SF pushes 
the particle off the geodesic, and thus causes orbital evolution which may be of practical importance. 

A second approach is based on the radiative Green's function [Q. Specifically, one can write the retarded Green's 
function G"'°' as the sum of two terms, namely 

where G^'^^ is the advanced Green's function. The first term is time symmetric, and consequently does not include 
the radiative part of the field. Instead, it relates to the non-radiative Coulomb piece of the field, which is the source 
for the divergence. If one considers then only the second term, namely the radiative Green's function, one can obtain 
an expression for the field of the charge, which is finite on the world line of the charge, and from this field obtain 
the SF. Although this approach is very successful in flat spacetime, it suffers from an inherent difficulty in curved 
spacetime. Specifically, in curved spacetime it is anti-causal: The radiative Green's function includes the advanced 
Green's function, which in curved spacetime has support inside the future light cone. Consequently, the momentary 
force on the particle depends in principle, according to the radiative Green's function, on the entire future history of 
the particle. When the future history is completely known, e.g., an eternal static particle, or eternal circular motion, 
this approach is expected to yield correct results. However, in general the future history is unknown, and might also 
be subjected to free will, such that an approach built on the radiative Green's function is unsatisfactory. An approach 
which is based solely on the causal retarded field is clearly preferable. (In what follows we shall indeed base our 
calculations on the retarded field.) Also this method ignores the conservative piece of the SF, which is included in 
the discarded time-symmetric part of G^^^. 

This approach was used by Gal'tsov to obtain the SF on scalar, electric, and gravitational charges in the spacetime 
of a Kerr black hole [Q. In particular, Gal'tsov found the SF acting on charges in a uniform circular orbit around a 
Kerr black hole in the weak field limit, and studied also the forces on static charges in the same limit. 
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In this Paper we focus on the SF acting on a scalar charge. Although a scalar charge is just a toy model for the more 
interesting and more realistic gravitational charge (a point mass), it already involves much of the properties of more 
realistic fields, especially the tail part of the SF. Yet, some of the complications associated with the gravitational SF 
are not invoked, notably the gauge problem of the SF. Also, the scalar field Green's function has just one component, 
which simplifies the analysis. Thus, despite its relative simplicity, the problem of the scalar field SF captures the 
essence of the physics, while avoiding some technical complications. For that reason the scalar field SF has been a 
very useful toy model. 

In what follows we shall compute the self force to linear order in the particle's field [i.e., to order (charge)^, and 
neglect corrections of order (charge)'* or higher. We thus assume that the particle's charge is much smaller than the 
typical length scale for the gravitational field, specifically the black hole's mass.] 

For a number of simple cases, e.g., static electric [|| or scalar |^ charges in the spacetime of a Schwarzschild black 
hole, or a static electric charge on the polar axis of a Kerr black hole the particle's field is known exactly in a 
closed form. Indeed, these exact solutions were used to find the SF for those cases However, in most cases 



an exact solution for the particle's field is unknown. For sufhciently simple spacetimes, e.g., that of a stationary and 
axisymmetric black hole, one can decompose the field into Fourier-harmonic modes, which can be obtained relatively 
easily. Specifically, the individual modes of the field (or the Green's function) satisfy an ordinary differential equation 
(whereas the field itself satisfies a partial differential equation), which for almost all cases can be solved (at least 
numerically using standard methods). In addition, it turns out that the individual modes of the field are continuous 
across the particle's world line and the resulting contributions to the SF of the individual modes are bounded. The 
divergency arises only at the step of summation over all modes. 

This prompted Ori to propose a calculation of the SF effects which is based on the retarded field and on a mode 
decomposition ||l^ (in that case the adiabatic, orbit integrated, evolution rate of the constants of motion in Kerr). 
More recently, Ori proposed to apply that method directly for the calculation of the SF . The greatest challenge, 
as was already mentioned, lies with an appropriate prescription for the regularization of the mode sum. When the 
individual modes of the SF are summed naively, the result typically diverges. The reason for this divergence is that 
the modes do not distinguish between the tail and the instantaneous parts of the SF, and contribute to both. [This 
occurs already in the case of a static scalar or electric charge in flat spacetime, when the position of the charge docs 
not coincide with the center of the coordinates. In that case the contribution to the SF of each mode (after summation 
over all azimuthal numbers m from — / to /) is independent of the mode number I, and is given by —q^/{2r'^), q being 
the (scalar or electric) charge, and r being the position of the charge. Obviously, the sum over modes diverges, and 
should be removed by a certain regularization prescription.] A mode sum regularization prescription (MSRP), which 
handles this divergence, and which is an application of the approaches of Mino et al [^ and of Quinn and Wald [|, 
was proposed by Ori p3| , p^ . 

Next, we describe the MSRP very succinctly. Further details are presented in Refs. [|l^,^. The contribution to the 
physical SF from the tail part of the Green's function can be decomposed into stationary Teukolsky modes, and then 
summed over the frequencies w and the azimuthal numbers to. The tail part of the SF equals then the limit e — > 0" 
of the sum over all I modes, of the difference between the force sourced by the entire world line (the bare force '^^"^/^) 
and the force sourced by the half-infinite world line to the future of e, where the particle has proper time r = 0, and 
T = e is an event along the past (r < 0) world line. Next, we seek a regularization function h''^ which is independent 
of e, such that the series J^iC'^^'^fli ^ converges. Once such a bmction is found, the regularized tail part of the 
SF is then given by 

/T + t CO 
dr'W^G^^'iz^ir), z"(t')] = ^ C^^^V/, - h^) - d, , (2) 
;=o 

where is a finite valued function. MSRP then shows, from a local integration of the Green's function, that the 
regularization function h^^ — Qf^l + bf^ + c^/^*. For several cases, which have already been studied, MSRP yields the 
values of the functions a^,5^,c^ and analytically. Alternatively, a^,6^, and (but not d^) can also be found 
from the large-^ behavior of '^''™/^. As X^l^o ^'^^^ fji ~ ^m) converges, it is clear that the large-Z behavior of ^'^'^'^fj^ is 
identical to the large-^ behavior of /ij,. For more details on MSRP, and in particular on the local integration of the 
Green's function and the analytical derivation of the MSRP parameters, see Refs. 

The MSRP has been applied successfully for a number of cases, including the SF on static scalar or electric charges 
in the spacetime of a Schwarzschild black hole [^ , a scalar charge in uniform circular orbit around a Schwarzschild 
black hole [^ , an electric charge in uniform circular motion in Minkowski spacetime , a scalar charge which is in 
radial free fall in the absence of the SF into a Schwarzschild black hole [|l^ , and the SF on static scalar or electric 
charges inside or outside thin spherical shells [^^. In all these cases the MSRP parameters a^,b^,c^ and d^ were 
known analytically, and were used in the regularization of the SF (for a list of these parameters see Ref. UM). The 
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SF has not been calculated using MSRP yet for cases where the MSRP parameters are unknown. However, there are 
interesting cases where these parameters have not been found analytically yet. Specifically, they have not been found 
analytically for any spacetime which is not spherically symmetric. In addition, when the spacetime of a spinning black 
hole is concerned, their analytical derivation is expected to be considerably more complicated than in the spherically 
symmetric case: the Green's function is generally time dependent, and for the corresponding wave evolution in the 
spacetime of a rotating black hole different I modes of the field couple |Q , a phenomenon which docs not happen in 
spherically-symmetric spacetimes. As was noted above, the MSRP parameters and c^, can be found also by 

studying the large-/ behavior of the individual modes of the bare force. However, without knowledge of the MSRP 
parameter any regularized result will not be unambiguous. It appears, then, that a local analysis of the retarded 
Green's function is unavoidable. Remarkably, it has been found that in all the cases for which the MSRP parameter 
is known, it satisfies a simple relation with the local part of the SF. Specifically, it has been shown that in those 
cases df^ equals the sum of the ALD force and the Ricci-curvature coupled piece of the SF, i.e., in the case of a scalar 
charge, 

d^L = {Ufi - Ufj.iiau") + ^'f'^t^ ■ (3) 

It was then conjectured, based on the particular cases for which it was found to be true, that Eq. (|^) is generally 
satisfied, at least for large classes of scalar charges ||l5|,^ . (In Ref. this conjecture applies only to static spherically- 
symmetric spacetimes. Here we expand the domain of validity of this conjecture to include at the least also the 
spacetime of a stationary, axially-symmetric black hole. We note that if this conjecture is found to be valid in general, 
it would not be unexpected.) If this conjecture is true in general, the full SF is given by just 
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1=0 



and all the terms appearing in Eq. can be found by studying the individual modes of the bare SF only. [Even if 
this conjecture is not true in general, Eq. (m still holds for the classes of cases for which the conjecture is true.] 

In this Paper we shall assume that Eq. (pFholds for the case of a static scalar charge in the spacetime of a stationary, 
axisymmetric black hole, and compute the SF for that case. By a static particle we mean a particle whose Boyer- 
Lindquist coordinates r, 6, and ip are fixed. Such a particle is static with respect to an undragged, static observer at 
infinity, but it rotates in the opposite direction to the black hole's spin with respect to a freely falling local observer. 
This is the first application of MSRP for a spacetime which is not spherically symmetric. In addition, we also use 
MSRP without prior knowledge of the regularization parameters. By comparing our result for the SF with known 
results for a Kerr spacetime in the weak field limit (which were obtained using independent methods in Ref. [Q), 
we shall, in fact, prove the validity of Eq. (||) in that limit. Moreover, by comparing the SF to the flux of angular 
momentum across the black hole's event horizon, we shall prove that the conjecture (||) is satisfied also for strong 
fields. 

The SF on static charges in the spacetime of black holes was considered in a number of works. The question 
to be asked then is the following: How is the external force which is needed in order to keep the particle static 
changed because of the self interaction of the particle? For the case of a static electric charge q in the spacetime of 
a Schwarzschild black hole of mass M, Smith and Will and Frolov and Zel'nikov found that there is a 

repelling inverse cubic force, which in the frame of a geodesic observer who is momentarily at rest at the position of 
the charge q is given by 

F.^q'^St, (5) 

r being the radial Schwarzschild coordinate. This force arises from the tail part of the full SF. The Schwarzschild 
spacetime is Ricci flat, such that the coupling to Ricci curvature does not contribute to the SF. Also, a static charge 
in Schwarzschild has zero ALD force, such that the fuU SF (|) is given by just the tail force. 

For the case of a static scalar charge (where the scalar field is minimally-coupled) in Schwarzschild, Wiseman ^ 
found the interesting result, that the SF is zero (this was found earlier by Zel'nikov and Frolov as a particular 
case in the more general study of non- minimally coupled scalar fields). Since the SF again is given by just the tail 
part (as Schwarzschild is Ricci flat, and the ALD force for a static scalar charge is again zero), it turns out that the 
tail part of the SF for a static scalar charge in Schwarzschild is zero. This vanishing result is in some sense surprising: 
spacetime is curved in a non-trivial way (which indeed leads to a non-zero tail force for an electric charge) , such that 
an exactly vanishing result for the tail force is not intuitively expected. The following question then arises: Is this 
result just a consequence of the particular symmetries imposed, which cause the SF to take a zero value? When these 



4 



symmetries are relaxed, does the zero result for the SF persist, or do we find a non-zero result? Some aspects of 
the symmetries were indeed relaxed in subsequent works. In Ref. |l7| the scalar charge was not considered anymore 
as static. Instead, a scalar charge in uniform circular motion around Schwarzschild was considered, and indeed a 
non-zero tail force, proportional to the angular velocity squared, was found. This tail force vanishes in the limit of 
zero angular velocity (a static particle). In Ref. [ pO| the horizon condition was relaxed, and instead the spacetime 
was taken to be that of a thin spherical shell. Again, a non-zero SF was found, which vanishes in the limit that the 
radius of the shell approaches its Schwarzschild radius. In these cases, too, the non-zero SF arises from the tail part 
of the SF. We thus see that the zero result for the SF (and in particular for the tail part of the SF) for the case of a 
static scalar charge in Schwarzschild is just a degenerate case when more complicated cases are considered, namely, a 
scalar charge in circular orbit, or a spacetime which is Schwarzschild, but not that of a black hole. Is the result of a 
zero SF for a static scalar charge outside a Schwarzschild black hole then just an isolated result, or is it a particular 
case of a wider class of spacetimes? 

We shall consider this question for a number of generalizations of the Schwarzschild spacetime. First, we shall add 
to the black hole electric charge, thus making it a Reissner-Nordstrom black hole. The Reissner-Nordstrom spacetime 
is electrovac, and hence Ricci curved. However, it turns out that the Ricci part of the SF for a static scalar charge 
vanishes. Also the ALD part of the SF vanishes for that case. What about the tail part of the SF? We find that the 
tail part of the SF is zero, and so is the full SF. We thus see that by adding electric charge to the black hole, we 
do not change the zero result for the SF. What about adding angular momentum to the black hole? We study then 
the SF on a static scalar charge in the spacetime of a Kerr black hole. Spacetime is Ricci fiat, but there is a local 
contribution to the SF from the ALD part of the SF. This contribution is in the dj dtp direction (in Boyer-Lindquist 
coordinates). (Incidentally, the force which is required in order to hold the particle fixed in the absence of the SF 
has components only in the 8/ dr and 8/ 80 directions.) When we consider the tail part of the SF, we find that it 
is still zero. Namely, the full SF is given by just the local, ALD part of the SF. We then consider the most general 
stationary, axially symmetric black hole, by adding both angular momentum and electric charge to the black hole, 
turning it into a Kerr-Newman black hole. Spacetime now is Ricci curved, and there is a non-zero contribution to 
the SF from the local Ricci-coupled part of the SF. There is also a contribution from the local ALD part of the SF. 
However, when we consider the full SF, we find that it equals just the sum of the two local terms, such that the tail 
part of the SF is again zero. 

We thus conclude that the tail part of the SF on a static scalar charge is zero for all stationary, axially symmetric 
black holes. The zero result in the spacetime of a Schwarzschild black hole turns out then to be just a particular case 
of a much wider class of spacetimes. In Schwarzschild also the local parts of the SF turn out to be zero, such that 
the full SF vanishes. However, in more general spacetimes the local terms are non-zero, such that there is a non-zero 
SF, but the interesting (and the difficult to find) part of the SF is the tail part, and it is zero for all stationary and 
axisymmetric black holes. 

We find that the full, regularized SF on a static scalar charge q in the spacetime of a Kerr-Newman black hole with 
mass M , spin parameter a, and electric charge Q, is given in Boyer-Lindquist coordinates by 

Here, q is the scalar charge of the particle, the horizon function A = — 2Mr + + Q^, and Yj = r"^ + a? cos^ 9. 
Equation (|^) for the full, regularized SF acting on a static scalar charge in the spacetime of a stationary, axisymmetric 
black hole is our main result in this Paper. 

One striking feature of the SF is that it diverges as the static limit (the outer boundary of the ergosphere) is 
approached. This situation is much different from that of a static electric charge in Schwarzschild, given by Eq. (|^), 
where the SF is bounded. Also, the ratio of the SF in that case to the regular force which is needed to be applied 
in order to hold the charge static (in the absence of the SF) tends to zero as the black hole's horizon is approached. 
Indeed, Smith and Will found that the SF is just a tiny correction, which becomes important only when the theory 
is no longer expected to be accurate, namely, when quantum effects are needed to be considered p^ ]. (Smith and 
Will found that, at its maximum at r = 3Af , the acceleration due to the SF becomes comparable with the regular 
acceleration if the Schwarzschild radius of the black hole is smaller than the classical electron radius, and the distance 
of the electron from the horizon is smaller by two orders of magnitude than the electron's Compton wave length.) 
In the strong field regime we thus find the SF acting on a scalar charge held static in the spacetime of a stationary, 
axisymmetric black hole to grow unboundedly as the static limit is approached. Note, that also the regular force 
which is needed to keep the particle static in the absence of the SF diverges in the same limit. The ratio of the 
acceleration due to the SF, a^^, and the regular acceleration, a™^^ for a Kerr black hole and on the equatorial plane, 
is found to be 
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(7) 



where q is the scalar charge of a particle with mass fj,, and M and a are the mass and spin of the black hole, 
correspondingly. By the accelerations of the left hand side we mean the corresponding magnitudes. As the static 
limit is approached (on the equatorial plane the static limit is at r = 2M), this ratio diverges, which signifies that 
the acceleration due to the SF becomes dominant over the regular acceleration. Also, there is a finite value of r for 
which the two accelerations become equal. When the charge-to-mass ratio q/n of the scalar particle is large, the two 
accelerations become comparable at large distances from the black hole. 

The organization of this paper is as follows. In Section |l^ we derive the equations governing the field, and obtain 
the expression for the SF. In Section III we derive the SF analytically in the weak field regime, and in Section IV we 
evaluate the SF numerically in the strong field regime. In Section Mwe derive the SF using the far field and balance 
arg uments, and compare our results with those obtained in Sections |III| and IV using the near field. Finally, in Section 
VI we discuss the properties of the SF. 



II. FORMULATION 



Consider a static scalar charge in the spacetime of a stationary, axially symmetric, black hole, i.e. the Kerr-Newman 
black hole. By a static particle we mean here that the particle's spatial position is fixed in Boyer-Lindquist coordinates. 
The background spacetime is described by the Kerr-Newman metric, which in Boyer-Lindquist coordinates assumes 
the form 



2Mr - Q2 



dt^ -^dtdip + —dr^ + T^de^ 

E A 



(8) 



where E = -|- cos^ 6*, A = + a'^ + - 2Mr, and C = {r^ + a'^)^ + a^{2Mr ~ Q^) sin^ 9 ^ (r^ + a^)^ - Aa^ sin^ 6*. 
We use units in which G = c = 1 throughout. Here M, a and Q are respectively the mass, angular momentum per 
unit mass and electric charge of the black hole. 

The linearized field equation for a minimally coupled, massless scalar field $ is given by 

V^V''$(a;") = -47rp(x") , (9) 

where denotes covariant differentiation compatible with the metric (^. The scalar charge density p is given by 

_6^x^' - z^(r)] 



dr- 
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(10) 



Here q is the particle's total scalar charge; t is the proper time; g = —T? sir? 9 is the metric determinant; and 
is the world line of the charge. We assume that the charge is placed at a position (ro, 6*0, (/?o)- To solve Eq. (||), we 
decompose p and into a sum over spherical harmonics Yi„i{9, ip): 
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(11) 
(12) 
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where Y,o — + a? cos^ ^o. Substituting Eqs. (pl|) and (|T2|) into Eq. (0), we find 



A0™ +2(r-Af)07-f 



A 
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2Mro - g2 



So 



Yhn{0o,'Po)Sir - ro) 



(13) 



where commas denote partial derivatives. The boundary conditions for $ are that $ vanishes as r — *■ cxd, and is 
regular on the future event horizon. Regularity of $ on the future horizon is equivalent to the field being derived from 
the retarded Green's function. (Regularity on the past event horizon is similarly related to the advanced field.) The 
solution of Eq. ([T^ ) can be written as 



6 



where 



2Mro - 
So 



(14) 



(15) 



Here Aq = rl+a'+Q''-2Mr^, r> = max(r,ro), r< = min(r,ro), 4™](ro) = </>i'"(ro)4";(ro)-4'"('^o)'/'i" ('^o) 

is the Wronskian determinant evaluated at r = tq; 0^™ and (/>2™ are two independent solutions of the homogeneous 
equation 



/^(t>^rr + 2(r - M)(/)_r + 



-/(/ + 1) 



= , 



(16) 



with (pl^"" satisfying the boundary condition at infinity, and (f)^^ chosen appropriately to make <i>'™ regular on the 
future event horizon. 

Equation (|l^) has three regular singular points at r_|_ , r_ , and at infinity, where r± = M ± ^ — — a'^ are 
the outer and inner horizons of the black hole. We next move the regular singular points of Eq. ( |l6|) to ±1,00. This 
is done by the transformation 



Equation (Rq) then becomes 



(1-z^ 



z(r) 



2r — r+ — r_ 



1(1 + 1)- 



, 



(1-Z2)_ 

which is the associated Legendre equation. Here the degree /i is purely imaginary and is given by 

a 



(17) 



(18) 



(19) 



Two linearly independent solutions are Pi^{z) and Qf{z), the associated Legendre function of the first and second 
kinds ||25| , respectively. The functions Pi {z) and Qi{z) can be expressed in terms of hypergeometric functions: 



Pt{z) = 



1 



r(i - m) 



z + l 



m/2 



1 - z 



The hypergeometric function has the Gauss series representation 

p . , N {a)n{b)n r 

,F,{a,b;c;z) = l^-^z 



|z|<l. 



where Pochhammer's symbol is defined to be 

(x)o = 1 and (a;)„ ~ x{x + 1) • • • (x + n — 1) {n > 1) . 



(20) 
(21) 

(22) 
(23) 



Note that the hypergeometric function which appears in the Pj^{z) expression ( po|) is just a polynomial of order I, 
and the series expansion (|2^) is valid even though the magnitude of the argument is greater than one. 

As r ^ oo (z ^ oo), Pj^{z) ~ and Qf(z) ^ In order for $ to vanish at infinity, we must have 

4>2^iz) = Qfi^) up to an arbitrary multiplicative factor [Note, from Eq. (p^), that the value of t/)'™ does not change 
if or (/jj™ is multiplied by a factor independent of r] . 

To determine (/)5™(z), we first consider the case to = 0. The function <i>"^(r, 6*, (^) — (j){r)Yio{9,ip) is independent of 
tp, such that the singularity of the coordinate (p on the horizon does not complicate the analysis. The general solution 
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of the homogeneous equation is the Unear combinations of the Legendre functions Pi{z) and Qi{z). On the horizon 
z = 1, Qi diverges but Pi remains finite, so we have (f>'i{z) = Pi{z) up to a multiphcative factor. When m 7^ 0, we 
find it more convenient to expand the general solution as a linear combination of Pl^(z) and ''(z) instead of Pj^{z) 



and Qi{z). Note that P^ ^{z) is a solution of (18) since the equation is invariant when /i is changed to — /_*. The two 
functions behave as 

1 / r 4- 1 \ 1 /^_^ ±ima/{r+-r-) 



as r ^ r+ . Both solutions oscillate near the horizon because of the coordinate singularity of ip there. To remove this 
coordinate singularity, we consider the coordinate (p{ip, r) defined such that 

d(fi — dip + -rdr . 



A 



Upon integration, we have 



a 



^ = ^+ I ^dr = ^ + — -— In(--i) (25) 



plus an integration constant which we set equal to zero without loss of generality (it only relates to the origin of the 
coordinate ip.) Hence 



In terms of the [r,6,ip) coordinates (the ingoing Kerr-Newman coordinates), the two functions Pl^{z)Yim{0,ip) and 
Pi^^{z)Yi,n{S,(p) become 

Yi„,{e, ^) ^— ^ and — — — 1^„(0, ^) 



^-'yr-r+J r(i + M) 

as r — > r+ {z 1), respectively. (When regularity of $ on the past event horizon is required, the appropriate 
coordinate transformation is given by dip = dip — ^dr.) The first expression still oscillates near the horizon while the 

second one is regular. Combining with the to = case, we conclude that (f>^{^{z) = P^'^i^z) up to a multiplicative 
factor. 

Note that each mode of the potential <i>''"(r, 6*, (ys) = (fJ-™- {r)Yim{d , ip) is complex in general. However, it is easy to 
show that (j)''"™ = (<i)'™)*, such that the TO-sum in Eq. ( |T^ ) is real and so is the scalar field <f>. 
The modes of the bare SF, are given by 

/m= E <l^^.'^'"\r^,0^,V^) ■ (26) 

m— — / 

It follows from Eqs. ( p^ ) and ( p^ that and <i>'™ are continuous at the position of the charge, but <i>'™ is not. 
Hence fl is not uniquely defined. However, the MSRP guarantees that the regularized SF does not depend on which 
derivative [$'™(r|^), $'™(rg") (where are the one sided limits of r ^ tq from above or below, correspondingly) or 
the average of the two] we use [ p5| . In practice, we use the average derivative, i.e., we define 

/'-^ E i[^':r{rt) + ^'"\r^)\ ■ (27) 



2 ^ 

m— — / 



Explicitly, we find that 



/^(ro, 00, ^0) = i E w J>f^I°iu A W{r^)4>'Z{r^) + '/'^.(ro)^^ (ro)] F,:„(0o, ^o)Fz™(0o, V'o) , (28) 
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qS{ro) 



and 



7n— — l 



Wr 



0^™](ro) 



(r-o)0r 



(29) 



(30) 



Typically, /^'^ diverges when summed naively. We then use MSRP to regularize the SF as outlined in Section I. 
Recall that the MSRP parameters a^,b^, and can be determined by two independent methods. Specifically, they 
can be found by either (i) study of the large-/ behavior of the individual modes of the bare SF, or by (ii) a local 
analysis of the Green function. In this paper, we shall study the regularization parameters using the former method. 
(It is hard to apply the latter method here because of the following reason. Any time-dependent evolution of the 
wave equation in the spacetime of a rotating black hole has to handle mode couplings, and the Green's function is 
obviously time dependent [p2|.) 



We shall carry out the regularization procedure analytically in Section III to study the SF in the weak field regime, 
and then numerically in Section IV in the strong field regime. 



III. COMPUTATION OF SELF-FORCE IN THE WEAK FIELD REGIME 



In this Section, we consider the case for which the scalar charge is far away from the black hole, i.e. tq ^ M 
[zq = (2ro — r+ — ?'_)/(r+ — r_) 3> 1]. We can thus expand in Eq. (|2^) in powers of r^^ . 
The large zq expansion of 0'/™ and (j)'^ can be carried out by using Eqs. (EOh and (pl|). Specifically, we choose 



■^1 (^) = 7^r(l + M + l)Pr{z) ^ , ^ „x M ^ 2F1 -/,/ + !; l + /i;^ (31) 



(20! ' ' ' ' ^ ' ^ ' (2;)! r(i + /i) - 1/ V ' ' 2 

^ri.) = j\) Qf (^) = - rr'\F, (1 + '-±1^, . ^ I . ,32) 



It can be shown, from Eq. (|22|), that the above expressions are equivalent to 



) 5^ n! (-2/)„ yT~z, 
(z) = J- r 1 - 1 V'T Sll^UZ^A. (34) 



<^'r(-) = ( 1 + -V "7i - E ^'?:!'o,r ( ) (33) 

ti/2 00 fl + 

Up to this point, no approximation has been made, but Eqs. (|lj), (|33|), ( |34[ ) and ( |26| ) give us a convenient way to 
expand in powers of . Hereafter, we evaluate all quantities at the position of the particle. To simplify the 
notation we shall assume that the scalar charge is placed at (r, 6, (p) and drop all the subscripts "0" . 

We first consider the expansions for fl and fg. We shall compute up to the order r^^ and fg up to the order 
. Straightforward calculations give 

i^qS{r) fX,{Q,0) , 2[l{l + l)Xi{0,9) + 3^^Xi{2,6)] 1 

J r 



21 + 1 (2l-l){2l + l){2l + 3) z2 (2/-3)(2/- l)(2Z-Hl)(2/-h3)(2/-|-5) 

[l{l + + I - 3)Xi{0, 6) + 2(3/2 + 3/ - 5)72^,(2, 6) + 5j^Xi{4, 9)]^ + O ( ^] } (35) 



fo- q^NP - a2 - Q^S{r) \ _ ^^^^^'2^^ 1)^(2/ + 3) ^ + (2/ -^3)f2r-^?)f2/TlK2/ + 3)^2/^+ 5) + ^ ^ ^ ' ^^^^ 



where 



Xi{p,6)^ J2 mPy;,^ {9, p)Yi„, (0, ^) (37) 

ni— — l 
I 

^i{p,e)^ Y mPYC^{e,p)deYUO.^) • (38) 
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In writing Eqs. ( |35| ) and (|3^), we have used the results that ^((0, 9) = and Xi{p, 9) = 0) = ifp is a positive 
odd integer. We leave the detailed calculations of the two fu nction s Xi{p, 9) and £,i{p, 9) to Appendix Substituting 
the expressions of Xi and S^i from Appendix]^ [Eqs. ( |A17| )-( [A22| )], we obtain 



Sir) 
8tt 



1 



2/(? + l) 



(2;- l)(2Z + 3) 
;2 



1 



37^ 



so? 9 



6;(? + l) 



z2 (2;-3)(2;- l)(2Z + 3)(2Z + 5)' 



3) + (3/2 + 3/ - 5)7^ sin^ 9 + ^7-* sin^ 6I + ^-^^(1 - + 2) si: 

2 8 



v/A/2 - a2 - g272 sin 61 cos 6* | ^^^-^ 



1 



2/(? + 1) 



- 1)(2; + 3) z (2; - 3)(2; - 1)(2; + 3)(2Z + 5) ■ 



For large values of I we find 



qS{r) 



1 



37' . 2 
— — sm^ ( 

2 



1 3 

72 + s 



1 + 37^ sin^ 



15 



-7 sm 



167r 



VAi"^ - a2 - g2 72 sin 6* cos 6* 



1 1 

^+2 



1 + t7^ sin^ I 
4 



O 



O 



(39) 
(40) 

(41) 
(42) 



Hence we find that the MSRP parameters a,. = = up to the order of r~'^ (i-e., any deviation of a,, or from 
zero is of order or higher) [Note that S{r) gives an extra factor r~2]^ ag = ce = up to the order of (i.e., any 
deviation of ag or cq from zero is of order or higher), and and hg are given from Eqs. ([41]) and (p2|): 



g^(0 
Stt 



1 + ^ ( l + ^sin^^) 4 + ^ ( 1 + 372 sin2 0+^74sin4 0) 4 + ^ ( 4 
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and 



g5(r) 
IGtt 



)2 72 sin COS I 



- + o ( 1+47' sin' ( 



4 + 4 



[For the case for which the particle is on the polar axis of the Kerr-Ncwman black hole we find that 



1 - M/r 



2r2 l-2M/r+(a2 + Q2)/r2 



-1/2 



(43) 



(44) 



(45) 



We have checked the latter expression numerically in the strong field regime, and found complete agreement.] Assum- 
ing the conjecture (^), the r-component of the regularized SF is then calculate d by subtracting Eq. ( ^3| ) from Eq. ( |39| ) 
and then summing over I. The ^-component is evaluated similarly from Eqs. (^^ and (^). The results are 



00 , 



1=0 



1=0 

where we have used the fact that for any integer k 

00 - ^ 

1=0 



2{l - fc) + 1 2(/ + fc) + 1 



= 0. 



(46) 
(47) 

(48) 



Hence we conclude that any non-zero orthonormal r and 9 components of the SF [Ff- and Fg) are of order or 
higher. In the next section, we present strong numerical evidence to suggest that F^ and Fg are actually zero wherever 
the location of the charge. 

The (/^-component of the SF can also be computed in the same way. However, we find a better method to do the 
calculation, which is described in detail in Appendix H. We find that 
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ri ^ qSjr) 

J te -.9; 



r2 ) X! 



(21 + 1) 



(49) 



where (21 + 1)!! = 1 • 3 • 5 • • • (2/ + 1) and 



Kt(z) = 1 - 



m/2 



iFi 1 



— - A • 



(50) 



that the leading term of /' is of order r 



-21-2 



It follows from Eq. 

keeping terms up to r^^ , only finite number of terms with I < (N - 



f 

J IP 



when I > {N — 2)/2, such that the MSRP parameters — ~ 



So when we make an asymptotic expansion by 
2)/2 contribute. In other words, up to the order 
0. In the next Section, we show 



numerically that decreases exponentially with increasing I, which also suggests that no regularization is needed to 
compute F^. Expanding Eq. ([49|) in powers of r~^, we find that 



!=0 



r 2 



^) + 2a2(l-3cos2 9)]^ 



01 ^ 



(51) 



The leading order term of this expansion agrees with the result given in Ref. for the case where (5 = 0. This 
agreement implies that we use here the correct parameter d^, and that the conjecture (^ holds for the case studied 
here. Note that this conclusion is valid only for the weak-field regime of a Kerr spacetime. In the next sections we 
bring evidence that this conjecture holds also for cases for which Q 7^ 0, and also in strong fields. 



IV. COMPUTATION OF THE SELF FORCE IN THE STRONG FIELD REGIME 



Next, we compute the SF in the strong-field regime. To do so, we use the general expressions (p8|)-(pO[) for the modes 
of the field, and compute them numerically. In order to find the regularized SF, we shall use MSRP. As discussed 
above, we determine the MSRP parameters a^,b^, and by studying the large-Z behavior of the individual modes 
of the bare SF. We shall, however, check our numerical results for the regularization parameters in the sphe rically 
symmetric limit, where the regularization parameters are known analytically p5 1 . Th is is done in subsection IV A, 
where we study the case of a Reissner-Nordstrom black hole. Then, in subsection IV B| we consider the case of a Kerr 
black hole, and in subsection pVC| we consider the case of a Kerr-Newman black hole. We shall also compare our 



results for the strong field with the weak-field approximation of Section III 



A. Reissner-Nordstrom 



Our main goal in this Paper is to study the SF on a static scalar charge in the spacetime of a rotating black 
hole. We shall apply our computation also for the spherically symmetric, electrically charge, Reissner-Nordstrom 
black hole for two reasons. First, as we already noted, the SF acting on a static scalar charge in the spacetime of a 
Reissner-Nordstrom black hole has not been calculated yet. Our results in this subsection are thus new. Second, it 
will allow us to check our numerical code for a case where the MSRP parameters are known analytically. (Our code 
does not assume spherical symmetry, thus its computation of the angular dependence of functions is non-trivial.) 

It is clear from symmetry considerations that all the azimuthal components vanish, such that we check below only 
the radial component of the SF. In Reissner-Nordstrom, it was found analytically in Ref. that for a static scalar 
charge 

bR.N_ 1' ^~M/r 

2r2 1 - 2M/r + Q^r^ ' ^ ' 

and = = cj?^. Also, = 0. Figure ^ displays our results for a Reissner-Nordstrom black hole of mass M = 1 
and electric charge Q = 0.8M, for a particle at rg = 4M. In Fig. |^(A) we present the modes of the covariant radial 
component of the bare SF, for a Reissner-Nordstrom spacetime. These modes appear to approach a constant 
value in the large-Z limit. In order to check whether this limiting value coincides with 6^^, we plot in Fig. |l|(B) the 
difference between and b^^ as a function of /. As for Reissner-Nordstrom it was shown that the MSRP parameter 
df-^ = 0, and also the ALD and Ricci-curvature parts of the SF vanish, the conjecture (||) holds. Hence this difference 
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is just the modes of the regularized full SF. We find this difference to scale like / ^ for large values of I. This result 
is in full accord with the analytical results for the MSRP parameters. We also compute the regularized full SF. In 
Fig. |^(C) we plot — Z]j=o(/r ~ b?^) ^ function of I. We find that behaves like /"^ in the large-Z limit, such 
that — * as Z — > oo. We infer that Fr = Fl^°° = for a static scalar charge in Reissner-Nordstrom. Similar 
results were obtained for other choices of the parameters. We infer that both the tail, non-local part of the SF, and 
the local two terms of the SF vanish separately for a static scalar charge outside a Reissner-Nordstrom black hole. 
This generalizes the known result of a zero SF on a static scalar charge in Schwarzschild to Reissner-Nordstrom. 




FIG. 1. Self force on a static scalar charge in Reissner-Nordstrom. Upper panel (A): fl as a function of I. Middle panel (B): 
The difference fl — b^^ as a function of I. Lower panel (C): F}. as a function of I. The data presented here correspond to the 
parameters: r = 4M, a = 0, Q = 0.8M, and M = 1. For the actual computation we took 6 = n/4. 

B. Kerr 

When the black hole is endowed with non-zero spin, we no longer have analytical results for the MSRP parameters 
O/i, bp, and c^j. As noted above, these parameters can be found, however, from the large-Z behavior of the modes of the 
bare SF. Nevertheless, the MSRP parameter cannot be found by studying the modes of the bare SF alone, and a 
local analysis of the Green's function is necessary in order to determine it. Recently, a conjecture about the MSRP 
parameter was formulated |l^,|2^. According to this conjecture, the MSRP parameter equals the sum of the 
two local terms in the full regularized SF, as given by Eq. (|^). In the Kerr spacetime, which is Ricci flat, equals 
then just the ALD force. In the following we shall use this conjecture. The ALD part of the SF is given by 



F: 



ALD 



3^ (A 



sm 



' sin^ 6')5/2Ei/' 



(53) 



local 



= F, 



ALD 



We thus conjecture that — F°. Also, because spacetime is Ricci flat, F^ ^ 

Figure || displays the behavior of the individual modes of the bare SF for a static scalar charge in the spacetime 
of a Kerr black hole (i.e., the scalar charge has fixed Boyer-Lindquist coordinates r,9 and (p.) We choose here the 
parameters r = 2.2M, 9 = 7r/4, a = 0.2M, Q = 0, and M = 1, but similar results were found also for other choices. 
We find that and fg approach constants as Z ^ oo, and that the difference between two consecutive modes, fl — fl^^ 



^Kcrr „Kcrr 



and fg — fg ^, behave like I ^ in the large-/ limit. This behavior implies that the MSRP parameters a^°", 
and c^°" vanish. (Non-zero a^'^" implies linear growth of the modes with the mode number I, and non-zero c^'^" 
implies that the difference between two consecutive modes behaves like l^^. Moreover, non-zero cjf'"''' threatens the 
applicability of MSRP, as it implies a divergent d^'^" llll-) Recall that we have defined in the averaged sense. The 
corresponding "one sided" values are not zero, in general. We also find that fj^ decays exponentially with I for large 



values of I, which suggests that a 



Kerr 







with our results in the weak-field expansions in Section III 



c^'^", and, in addition, also b^'^" 



0. The vanishing of b^'^" is in accord 
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FIG. 2. Behavior of the individual modes of the self force for a static scalar charge in the spacetime of a Kerr black hole. 
Panels (A), (C), and (E): The r, 8, and (f covariant components of the individual modes of the self force, respectively (i.e., 
fr, fg, and /^), as functions of the mode number I. Panels (B), (D), and (F): The r, d, and covariant components of the 
difference between two consecutive modes of the self force, respectively (i.e., — fr^^ , fe — fg~^, and — f!p~^), as functions 
of the mode number I. The parameters for the data presented here are: r = 2.2M, 6 — n/'i, a = Q.2M, Q = 0, and M = 1. 



In order to compute the regularized SF we have to confront the difficulty of not having exact expressions for the 
MSRP parameter 6^^''" . However, is nothing but the hmit as Z ^ cxd of /^. We can thus approximate by 
studying the large-Z behavior of and extrapolate to infinite value of the mode number (e.g., through Richardson 
extrapolation). In practice, we can approximate b^°^^ by simply taking the mode of the bare SF with a mode number 
L much larger than the mode I up to which we sum over the modes to obtain the full, summed-over-modes, SF. That 
is, we compute the regularized SF according to 



CO 

E 

/ 



if. 



(54) 



(55) 



The functions 7?.'+^ and Sj^^ are defined by Eqs. and (|5^), respectively. Now, by definition 



j=i+i 



E 

3=1 + 1 



o(j- 



(56) 



where a;^ is the coefficient of the I ^ term in the I expansion of /' . In evaluating T?.'"*"^ we shall drop the 0{j '^) 



term in Eq. (p6|). This introduces an error which reduces like as I grows. By definition, the function £^ = 

'l2]=o(-ffi^ ^ /^f )' ^-iid inside the sum we just have a ^-independent expression, we find that £^ = (/ + 1)(/^ — f^)- 
Taking L 3> / S> 1, we approximate i^^j in practice by 



F. 



I 

E 



if. 



f ^ , 



(57) 



The error associated with this approximation has two contributions, i.e., 7?.'+^ and We evaluate 7?.'+^ 



1) Xfj^/l, where tp'-^^z) is the trigamma function. We also evaluated £^ « Ix^/L"^. Note, that can be evaluated 
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from the difference between two consecutive modes, i.e., a;^ « ^^(//i ^ ~ (when I is large enough). We find that 

£^/7?,Jj+^ ~ (V-^)^' such that when L ^ I the overall error is dominated by 7?,|j^^. 

Figure shows the behavior of the regularized SF. We choose here I — 330 and L = 1000. In Fig. ||(A) and |^(B) 
we show fl — br and as functions of I, respectively, and in Fig. ^(C) and ||(D) wc show fg — be and Fg as functions 
of I, respectively. As expected, we find both — br [Fig. ||(A)] and fg — bg [Fig. ||(C)] to behave like for large 
values of I. We also find that F^ [Fig. ||(B)] and Fg [Fig. g(D)] scale like for large values of I. When extrapolated, 
we conclude that as I ^ oo, F^ Fr ^ and Fg Fg — 0. 




u, -1 



log. 



1.5 
„ / 




log. 



1.5 

„ ; 



log. 
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FIG. 3. The regularized SF for a static scalar charge in the spacetime of a Kerr black hole. Panels (A) and (C): The r and 9 
covariant components of the individual modes of the regularized self force, respectively, (i.e., — br and fg — be) as functions 
of the mode number I. Panels (B) and (D); Fr = X/' =o(-^'^ ~ ^'') ^^'^ -^e ~ X/'=o(-^e ^ respectively, as functions of /. Panel 



(E) : The difference between i^^ = Ej=o(/^ " ^f) ^o^al SF, which is given by Eq. (^^, as a function of I. Panel 

(F) : as a function of I. The (unknown) values of br and bg were approximated by their respective values at L = 1000, i.e., 
by ff^^^^^ and /e^^^""", correspondingly. The parameters for the data presented here are: r = 2.2M, 9 — -k /A, a = Q.2M, 



= 0, and M = 1. 



We can approximate F^ and F^ even better by including approximate values for Tif^^^ and £f^. This is done in 
Fig. ||. In Fig. |^(A) we show F^ twice: without the inclusion of Tif^^ and [same as in Fig. 0(B)] and with their 
inclusion. We find that for large values of /, the latter behaves like l^^ (the former scales like I ). Similarly, in Fig. 
^(B) we show twice: without the inclusion of T?.^^^ and £g [same as in Fig. ||(D)] and with their inclusion. Again, 
we find that for large values of I, the latter behaves like (the former scales like Recall that we are using here 
only approximated values for 7?.|^+^ and £j^. For this reason, we do not expect their inclusion to yield an exact zero 
result for the r and 6 components of the SF. However, they do eliminate the leading order term in F^. and Fg, such 
that instead of a leading behavior we find an behavior. When this behavior is extrapolated to I — > oo, we 
again infer that F^ ^ Fr = and F^ ^ Fg = 0. 
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I°9l0 



FIG. 4. The regularized SF for a static scalar charge in the spacetime of Kerr. Upper panel (A): without the inclusion of 
711+^ and Sf- [same as in Fig. §(B)J (solid line) and with their inclusion (dotted line). Lower panel (B): Fg without the inclusion 
of 7^^+^ and £^ [same as in Fig. H(D)] (solid line) and with their inclusion (dotted line). 

In Fig. ||(F) we show as a function of I, and in Fig. ||(E) we show the difference between F^^ and F^'^'^^, which is 
given by Eq. (|53|). We find that F^^ approaches a non-zero value as Z — > 00. We also find that the difference between 
F^ and F^^'^^ decays exponentially in I for large-Z values. We infer that as I 00, F^ ^ F^ = F^'^'^^. We find similar 
results also for other choices of the parameters. 

When our results for the SF are combined, we find that the tail part of the regularized SF vanishes. That is, the 
tail part is given by 

drV.G*-' = ^(/^ - b,) - d, . (58) 

According to the conjecture (|^), = We find numerically that J^'jLoif^ ~ ^m) = Fj^'^^^- Consequently, we 

infer that F*^'' = 0. 



C. Kerr-Newman 



Next, we endow the black hole with an electric charge Q in addition to its spin parameter a. In this case, too, we do 
not have the analytical results for the MSRP parameters a^, 6^ and c^. As in the Kerr case above, these parameters 
can be obtained from the large-Z behavior of the modes of the bare SF. Again, we shall use the conjecture (^) regarding 
the MSRP parameter dfj,, according to which equals the sum of the two local terms in the full, regularized SF. 
As the Kerr-Newman spacetime has non-vanishing Ricci curvature, both terms contribute. We find that for a static 
scalar charge outside a Kerr-Newman black hole 



and 



Fa"" - -i^^aAsin^ 6 6^ , (60) 

3^ (A-a2sin2 0)3/2i]3/2 m' V y 



such that the total local piece of the SF is given by 

1 2 A • ■ 

3^ (A - a2 sin2 6')5/2Si/2 



1 Af2 /Q2 
rfocal ^ ^ 2^^ gij^2 g ^ /gl) 
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First, we study the behavior of the individual modes of the bare SF, which we compute numerically from Eqs. 
(^)"(^. Figure ^ shows the individual modes, and the difference between two consecutive modes, — fjj^^, 
as functions of the mode number I for /j, ~ r,9, f. In Fig. ^ we present the first 336 modes (i.e., ^ = 0, . . . , 335) for 
the parameters r = 2.2M, = 7r/4, a = O.IM, Q = O.IM, and M = 1, but similar behavior was found also for other 
choices for the values of the parameters. We find that the r and components of the modes approach a non-zero 
limiting value in the large-/ limit, as the difference between two consecutive modes behaves like for large values 
of I. This behavior of the individual modes implies that the MSRP parameters a^, a^, c^, and ce vanish. (Recall that 
we are using here the "averaged" value for and ae- The "one-sided" values are, in general, non-zero.) (Note, that a 
non-zero value for Cr or ce implies that the difference between two consecutive modes should scale like .) However, 
the MSRP parameters br and bg are non-zero. These parameters correspond to the I oo limit of the individual 
modes fl and fg, respectively (recall that and ae vanish). The (p component of the individual modes of the SF 
drops off exponentially for large values of I, and we infer that a^p — 0, b^p — 0, and = 0. This is in agreement with 



the weak-field expansion of Section III 





FIG. 5. Behavior of the individual modes of the self force for a static scalar charge in the spacetime of a Kerr-Newman black 
hole. Panels (A), (C), and (E): The r, 6, and ip covariant components of the individual modes of the self force, respectively 
(i.e., fl, fg, and /^), as functions of the mode number 1. Panels (B), (D), and (F): The r, 6, and cp covariant components of the 
difference between two consecutive modes of the self force, respectively (i.e., fi — fr~^, fl — /s^^i and — f!p~^), as functions 



of the mode number i 



The parameters for the data presented here are: r = 2.2M, 6 = 7r/4, a = O.lAf, Q — O.IM, and M = 1. 



Next we compare our results for the modes here, with our asymptotic expansions in Section III . Figure ^ displays 
the relative difference of the asymptotic values for the modes of the bare SF [given by Eq. (|3^) for the r component 
and by Eq. (Eo|) for the 6 component] and the full values of the modes of the bare SF [which we compute numerically 
from Eqs. (|2^) and (p9|)] as functions of r (in units of tsl = M + V — cos^ 9, the r value at the static limit) 
for different values of the mode number /. At large distances we find both the r and 9 components to agree with the 
asymptotic expansion at the right orders. Near the static limit, they of course disagree. Note, however, that even 
very close to the static limit the r component of the modes of the asymptotic expansion does not deviate from their 
full-expression counterparts by much more than 10%. We find similar results also for other choices of the parameters. 
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log,o(^/^SL) l°gio(^^^SL) 
FIG. 6. Comparison of the modes of the bare force as given by the asymptotic expansions [Eq. (|39|) for the r component 
and Eq. ( |4o| ) for the 6 component] and the full expressions which we compute numerically from Eqs. (^8|) and (^), respectively. 
We show the relative difference of the two expressions as a function of r/rsh for a number of mode numbers I. Left panel: the 
r component. Right panel: the 9 component. In both panels we show the modes: I = 1 (solid line), 1 = 2 (dash-dotted line), 
/ = 3 (dashed line), and 1 = 4 (dotted line). The data here are shown for the following values for the parameters: Af = 1, 
a = 0.6M, Q = 0AM, and 6 = 7r/4. 

Figure ^ displays the behavior of the regularized SF. We choose here again I = 330 and L = 1000. In Fig. 0(A) 
and 0(B) we show — br and as functions of /, respectively, and in In Fig. 0(C) and 0(D) we show fg — be and 
Fg as functions of I, respectively. Similar to the Kerr case, we find both /; - br [Fig. 0(A)] and fg - be [Fig. 0(C)] to 
behave like for large values of I. We also find that F^ [Fig. 0(B)] and Fg [Fig. 0(D)] scale Hke l^^ for large values 
of I. When extrapolated, we conclude that a,s I oo, F^ Fr — and Fg Fg — 0. 
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FIG. 7. The regularized SF for a static scalar charge in the spacetime of a Kerr-Newman black hole. Panels (A) and (C): 
The r and 6 covariant components of the individual modes of the regularized self force, respectively, (i.e., fl — for and /g — bg) as 
functions of the mode number I. Panels (B) and (D): F^i = X/j=o(-^'^ ~ ^'") ^^'^ ^t> ~ 'Yl!j=o^fe ^ respectively, as functions 
of /. Panel (E): The difference between = Y.]=o^f^ " ^-fi) ^-nd the local SF, F^"''''', as a function of I. Panel (F): F^ as a 
function of 1. The (unknown) values of for and bo were approximated by their respective values at L — 1000, i.e., by Z^^^""" 
and /e^^^""", correspondingly. The parameters for the data presented here are: r = 2.2M, 6 = tt/4, a = O.lAf , Q = O.IM, and 
M = 1. 



In a similar way to our analysis of the case of a Kerr spacetime, we can approximate and Fg even better by 
including TZ''^^ and f^. This is done in Fig. ||. In Fig. ^(A) we show F^ twice: without the inclusion of TZ^:^^ and £^ 
[same as in Fig. 0(B)] and with their inclusion. We again find that for large values of I, the latter behaves like 
(the former scales like Similarly, in Fig. H(B) we show Fg twice: without the inclusion of TZ'g^^ and £g [same as 

in Fig. 0(D)] and with their inclusion. Again, we find that for large values of I, the latter behaves like (the former 
scales like Recall that we are using here only approximated values for TZ''^'^ and For this reason, we do not 

expect their inclusion to yield an exact zero result for the r and 6 components of the SF. However, like in the Kerr 
case, they do eliminate the leading order term in F^ and Fg, such that instead of a leading behavior we find an 
l~ behavior. When this behavior is extrapolated to / —^ oo, we again infer that F^ ^ Fr = and F^ Fg ^ 0. 




_gi , 1 , , , , , , 1 
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FIG. 8. The regularized SF for a static scalar charge in the spacetime of a Kerr-Newman black hole. Upper panel (A): F'. 
without the inclusion of TZ'^^ and £^ [same as in Fig. 0(B)] (solid line) and with their inclusion (dotted line). Lower panel 
(B): Fg without the inclusion of TZ^g^^ and £g [same as in Fig. M(D)] (solid line) and with their inclusion (dotted line). 

In Fig. 0(F) we show F^ as a function of and in Fig. 0(E) we show the difference between and F^^'^^, which is 
given by Eq. (|6l|). We find that approaches a non-zero value as Z ^ oo. We also find that the difference between 
F^ and F^Jf'^^' decays exponentially in I for large-Z values. We infer that as Z ^ oo, F^ ^ F^ = F^'^^^. We find similar 
results also for other choices of the parameters. 

As in the case of a Kerr spacetime, when our results for the SF are combined, we find that the tail part of the 
regularized SF vanishes. 



V. FAR FIELD COMPUTATION OF THE SELF FORCE 



In Sections HI and IV we computed the SF by using the field (and its gradient) evaluated on the particle's world 
line, i.e., by using the near field. In this Section, we shall compute the SF using the far field (evaluated asymptotically 
at infinity and at the black hole's event horizon), and demonstrate the compatibility of the two approaches. For 
simplicity, we shall restrict our considerations in this Section to the case of a Kerr black hole. Specifically, we shall 
show that the covariant t and (p components of the SF {Ft and F^p, respectively) in the case of a Kerr black hole 
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can also be inferred from balance arguments pertaining to the global conservation of energy and angular momentum. 
Specifically, we shall deduce Ft and by calculating the fluxes of energy and angular momentum, associated with 
the charge's scalar field, flowing out to infinity and down the black hole's event horizon. We shall show that the 
results of these far-field calculations agree with the near-field calculations we have performed in Sections III and 
In fact, by showing this agreement we demonstrate the applicability of the MSRP, and the validity of the conjecture 
to the problem of interest. 

Since we need to evaluate the scalar field near the horizon and at infinity, we place the particle at (ro,Oo,ipo) in 
this Section to avoid confusion. 

The rate of change of the particle's four- momentum due to the SF, F^, is given by 



dT 



(62) 



where r is the proper time. In particular, we have 



1 dL 



v^-gtt dt 



(63) 



where L = is the angular momentum of the particle along the black hole's rotation axis, and we have used the 
fact that dt/dr — 1/ ^—gu (recall that the charge is static). The sum of the rate of change of the particle's angular 
momentum (i.e., dL/ dt), the total amount of angular momentum (per unit time) flowing out to infinity {J^^) and 
down the black hole's event horizon (.^^jfoio) must be zero, if global angular momentum were to be conserved. Hence, 



- dL/ dt = T^+ 



hole 



(64) 



The value of T]^ is given by (see Chapter 5 of Ref. 123]) 



= lim / TJ & N dfl 



(65) 



where ^i-^) = d/dcp is the axial Killing vector, dfl = sinO d9 dip, and the stress-energy tensor T^j, associated with the 
scalar field $ is given by 



1 

47r 



(66) 



The total angular momentum flowing down the event horizon per unit time, ^i^^ic, is given by |2£ 



J^t.o = ^Hm J [~2Mt+T^, Ihh] , 



(67) 



where is one of the basis vectors of the Hawking-Hartle tetrad. It is an outgoing tetrad which is made well 



behaved on the future event horizon [pl| , p0y29[| . The components of I'^^j, in Boyer-Lindquist coordinates {t,r,9,ip), 
are given by 



''HH 



,0, 



2'2(r2 + a2)' '2(r2+a2) 



(68) 



Note that although some of the components of T^^, in Boyer-Lindquist coordinates, diverge on the event horizon, 
•^hoie remains finite. In fact, the divergence is due entirely to the coordinate singularity. However, J^^^i^, being a 
scalar, is independent of the choice of coordinates. 

The asymptotic e xpr essions of the scalar field <I> at infinity and on the horizon can be deduced from Eqs. (p^), (pH), 
(H), (|3), ( ^TTl) , ( |B8|) , (|l|), (|l|) and (H). The results are 



<S>{r,d,ip) 
<i>ir,e,ip) 



1 - 



2Mro q 



Yim{0,(p) r ^ r+ 



(69) 
(70) 



where is defined by dr^,/dr = (r^ + a2)/A, km = —mujj^ = —ma/(2Mr^), and Zim is given by 
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Zlm — dim 



VA'P - S{ro) (1 - z^) 



z'+^ (2/ + 1)!! 



+ im'-f) 



(71) 



Here C;™ are complex constants ol unit modulus, i.e. |C;„i| = 1. Substituting Eq. (|69| ) into Eqs. (|66| ) and (|65|), we 
find that = 0. This is expected since the particle is static relative to static observers at infinity, so no radiation 
is emitted. The rate of change of the particle's angular momentum is then given by 



dL 
~dt 



r— >r+ 



(72) 



Straightforward calculations yield 



hole 



EE 

oo / 

EE 



in z, 



2/-2 



^0 - 



n(mv+j^) 



i^r(^o) 



(2; + l)!! 



(73) 



Yi:j6o,^o)YUeo,^o) ■ (74) 



The covariant (^-component of the SF, F^, is then calculated by Eq. (|63|). Using Eqs. (jl^) and ([l7|), we finally obtain 



00 / 



EE 



qaS{ro) 



1 - ^ I 

^5, 



n(™v+/) 



(2^ + 1)!! 



(75) 



which agrees with Eq. ( fl^ ) (Recall that the MSRP parameters = = = 0). 

Similarly, we can deduce the covariant time component of the SF, Ft, by the energy balance argument. We have 



Ft = 



dpt 
dT 



1 dE 



V^gTt dt 



(76) 



where E = —pt is the energy of the particle. The sum of the rate of change of the particle's energy (i.e., dE/ dt), the 
total amount of the energy (per unit time) associated with the scalar field flowing out to infinity {J-^) and down the 
horizon {^^^i^) vanishes because of global conservation of energy, i.e.. 



dE/dt = T^+Tt,, . 



The relevant formulae are [p8'30|: 



■> hole — 



lim 

r — *oo 



^hm y 2Mr+ T^, ^["^^ d^ , 



(77) 

(78) 
(79) 



where = d/dt is a Killing vector. Straightforward calculations yield JF^ = = .?^i^ic- Hence dE/ dt = = Ft, as 
expected. 



VI. PROPERTIES OF THE SELF FORCE 

From our numerical study we infer that the full, regularized SF on a static scalar charge in Kerr-Newman is given 

by 

FfF = -q^aA sin2 9 -Q^ (80) 

" 3^ (A - a2 sin^ 0)5/2i]i/2 ^ ' 

We note the following properties of this result. 



20 



The SF vanishes as a ^ 0, as expected. In that Hmit the black hole becomes Reissner -Nord strom (or 



Schwarzschild in the lack of electric charge), for which the SF vanishes, as we found in Section [V A 



When t ^ —t the SF reverses its sign. [Notice that the SF has the form of a x (even powers of a).] This is indeed 
expected, because under time reversal the black hole reverses its spin, and rotates in the opposite direction. 
This change of sign under time reversal implies that this SF is dissipative. 



The SF diverges as A — a sin 9^0. This is not surprising, as A — a sin 9 = defines the static limit, beyond 
which (inside the ergosphere) no timelike static trajectories exist. Consequently, our problem of finding the SF 
on a static particle becomes ill posed beyond the static limit. 

The SF is only in the d/ dip direction. The addition to the regular force which is needed to keep the particle 
in its static position is, therefore, orthogonal to the regular force (which has components only in the d/ dr and 
d/ d9 directions). 

Denoting the magnitude of the acceleration of the static particle in the absence of the SF by a''^^, and the 
magnitude of the acceleration due only to the SF by a^^, we find the ratio of the two accelerations to be given 

by 



1 



asmt 



3 /i 



E1/2A1/2 



I A [M (r^ - cos^ 9) - rQ^] ' + ^a^ (2Mr - 



''29 



-1/2 



(81) 



where is the mass of the scalar charge. Clearly, this ratio diverges as the static limit is approached. In this 
sense, the SF is not a tiny correction for the regular external force which is exerted in order to keep the particle 
fixed: the SF becomes dominant over the regular acceleration when the particle is sufficiently close to the static 
limit. The origin of the divergence of this ratio is in the "damping part" of the ALD force. The "Schott part" 
of the ALD force diverges too approaching the static limit, but at a slower rate than the "damping part". The 
part of the SF which couples to Ricci curvature diverges at the same rate as a''°s approaching the static limit. 
The leading order divergence in Eq. (^Tj) comes then from the "damping part" of the ALD force. This can be 
readily understood by the following argument. The "damping part" of the ALD force diverges like a'^°^'^ , such 



that its ratio to a'''^ is expected to diverge like a'^°^. Indeed, 
static limit. 



diverges like (A - 



- sin^ 



9) approaching the 



The scalar field theory does not restrict the charge-to-mass ratio of a scalar particle, q/fi. (For other field 
theories we find the charge-to-mass ratios to span many orders of magnitude: it is unity for a gravitational 
charge, but 2 x 10^^ for an electron.) We can thus view q/fi as a free parameter in Eq. (^ij). Recall, however, 
that q/M is assumed to be a small quantity. As r ^ oo the ratio a^^/a'^°^ 0. Thus, there is a value of 



(outside the static limit) where a' 



SF 



.rcg 



(Recall, however, that the SF acceleration is in the 8/ dip direction. 



while the regular acceleration is in the 8/ dr and 3/ 89 directions. 
Kerr, we find that a^^ equals a"^^ at 



Specializing now to the equatorial plane in 



-M 



2^ 1 , ,1 
3 ^ 



27a- 



M 



16/iAf + 3W3a^ ( 27a 



M 



32/iAf 



27a— + 16fiM -h 3 W 3a-^ ( 27a— + 32//Af^ 
M y M \ M J 



(82) 



For small values of ag^/(^M^), this value of r can be expanded in powers of ag^/(/iM^). We find that 



288 (J j {Jlfi-S^'*" 



M 



M 



such that for small aq'^/{iiA'P) the two accelerations become comparable only very close to the static limit. 
Figure ^ displays this value for r vs. the free parameter q/fj, for three values of q/M. Keeping a/M and q/M 
fixed, we find that at small values oi q/ ji the two accelerations become comparable only very close to the static 
limit. For large values of q/ fx they become equal at values of r which scale like {q/nY^^. The change in the 
behavior occurs near q/ ^ ~ M'^/{aq). Recalling that q/M ^ 1 and that a/M < 1, we find that at the change 
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of behavior in Fig. ^, g//.t ^ 1. In order to have r ^ M (the distance from the black hole at which the two 
accelerations become comparable is very large) we should thus require qj [i"^ M"^ /{aq). 
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FIG. 9. The value of r/rsi, where a^^ equals a"^^ as a function of on the equatorial plane of a Kerr black hole for 
a/M = 0.5, for three values of q/M: q/M = 10"* (solid line), q/M = 10"^ (dashed line), and q/M = 10"^° (dotted line). 

• At large distances this SF agrees with the SF found by Gal'tsov ||^, i.e., at large distances we find that 

= ^asin^ (l + 3— + ^[3(51/2 - Q^) + 20^(1 - 3cos2 0)]1 + O (r-^)! , (83) 

3 \^ r 2 j,^ \ / j 

whose leading term agrees with Gal'tsov's result when Q = Q (a Kerr black hole). Notice, that this expansion 
coincides with the asymptotic solution we found above in Eq. (|5l|). 

• The direction of the SF is in the direction of the spin of the black hole. Namely, in order to hold the particle 
static, the apphed external force should be in the direction opposite to the spin. As noted by Gal'tsov 0], this 
direction can be explained as a tidal friction effect |^^: As the particle is accelerated in the direction of the 
black hole's spin, global conservation of angular momentum implies that the black hole is accelerated in the 
direction opposite to the spin, such that the black hole tends to spin down, and its rotational energy is being 
dissipated. 

• When the particle's position is off the black hole's polar axis, the black hole is immersed in an external field, such 
that the entire configuration is not axially symmetric. From Hawking's theorem |26| , stating that a stationary 
black hole must be either static or axisymmetric, it then follows that the black hole cannot remain stationary: 
it must evolve in time until it has become static or until it has achieved an axisymmetric orientation with the 
external field As the scalar field is stationary, there can't be any flux of energy down the event horizon, 
as measured by a static observer at infinity. Hence the black hole's mass M is unchanged. As the black hole's 
surface area A, given by ^ = %nM{M + ^/ M'^ — a^) must increase, it follows that its angular momentum a 
must decrease, or that the black hole spins down. The dissipated rotational energy of the black hole does not 
escape to infinity, as the field is strictly static there. Instead, it fiows down the hole as seen by a local, dragged 
observer. Specifically, a local observer who follows a timelike orbit must be dragged inside the ergosphere. Any 
such observer will see the field as time dependent, and will see a flux of energy down the event horizon, whose 
origin is in the black hole's rotational energy. However, when the particle is on the black hole's polar axis, the 
SF vanishes according to Eq. (pO|). In this case, the black hole is immersed in an axisymmetric field, such that 
it can remain stationary. The flux of angular momentum as viewed by a static distant observer vanishes, and 
the black hole's spin is unchanged. 

• When Newton's constant G and the speed of light c are re-introduced, and the SF (^) is expanded in powers 
of G/c^, wc find that 
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a sm 



l + 2(l-3cos^6') ^ 



3+ (5 - Mcos^ 9) ^ 



ays 



O 



G 



(84) 



such that we find that this SF is a (post) ^-Newtonian effect. Note, that at the (post)^-Newtonian order the 
effect has contributions both from the and the terms. We also write explicitly in Eq. ( ^4|) the (post)^- 
Newtonian order of the SF. 
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APPENDIX A: 

EVALUATION OF Xi{p,e) AND $,iip,9) 

When calculating the SF in the weak field region, we encounter the functions Xi{p, 9) and £,i{p, 9) defined as 

I 



Xi{p,9)^ mPY:„,i9,^)Yi„,i9,ip) 



m— — l 
I 



^i{p,9)^ mPYC^{9,^)-Yi^{9,^) . 



rn— — l 



(Al) 
(A2) 



The spherical harmonics are of the form Yim{9,Lp) = 3^;m(^')e*'"'^, where yim{S) is a real-valued function. It follows 
that Xi{p, 9) and 5i(p, 9) are independent of ip. So we can take = and write 



Xi{p,9)^ Y mP[YUO,0)f 

m—~l 
I 

^i{p,9) = Y mPYU0,0)-^Yi^i9,0) = ~Xi{p,9) 



(A3) 
(A4) 



m— — l 



The spherical harmonics have the property that Yi,^m{9,(p) = {~l)™Yi^{9, ip), which implies that [Yi^^rn{9,0)]'^ 
[Yim{d, 0)f. So for p > 1, we have 



XKp,0)= ^mP[l + (-l)^][l,™(0,O)] = 



(A5) 



It follows from Eqs. ( |AE| ) and ( [A4| ) that Xi{p, 9) = ^i(p, 6*) = if p is an odd integer. 
To evaluate X/, we use the addition formula for spherical harmonics: 



m— — L 



where Pi is the Legendre polynomial and 

1/ = cos 9 cos 9' + sin 9 sin 9' cos{(p — ip') . 



(A6) 



(A7) 



It follows from Eqs. (|A6|), (^ and P;(l) = 1 that ^((0, 9) = {21 + l)/(47r) and £_i{0, 9) = 0. To compute Xi for p > 1, 
we set 9' = 9, ip' = and differentiate Eq. (AC) p times with respect to (p. After some rearrangement, we obtain 
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Xi{p,9) 



,21 + 1 dP , 



in dtpP 



(A8) 



with v = cos^ 6 + sin^ 9 cos ip. For simpli city , we only compute Xi{p, 6) for p = 2 and p = 4 here. Generalization to 
other values of p is straightforward. Eq. (A8) gives 



Xi{2,e) 



^^sin2 0P/(l) 
47r ' ^ ' 

2^ + 1 , o 



47r 



sin" 0F/(l) + 3sin* 



To compute and Pl'{l), we use the recurrence relation 

= a;i^;(a;) + (n + l)P„(a:) . 

Differentiating the above equation, we have 

P:+,{x)^xP:{x) + {n + 2)K{x) . 

Evaluating the two equations at x = 1 and using the fact that -P,i(l) = 1, we obtain 

PUii^)-P:,{l)=n+l 
/^Vi(l)-^'(1) = (^^ + 2)P/(1). 



(A9) 
(AlO) 

(All) 
(A12) 



(A13) 
(A14) 



Since Po{x) = 1, we have Pq{1) = -Po'(l) = The difference equations (A13) and (A14) are then solved by summing 
over n on both sides from n = to (/ — 1). The results are 



2 

{I - + + 2) 



Combining our results, we finally have 



XK2,^)^^i(^sin^^ 



Att 2 
21 + 1 



^SttH sin2 e + -{l~ l)l{l + 1){1 + 2) sin^ ( 



The values of ^; (p, 9) can then be computed by Eq. ( A_4 ) . The results are 

6(0,^) -0 

, , 21 + 1 1(1 + 1) 
(i{2,0) = —— ^ : ' singcosf 



6(4,0) 



47r 2 
21 + 1 



■ smc/cosf 



fc^ + J(/-l)/a + l)(/ + 2)sin2 



(A15) 
(A16) 

(A17) 
(A18) 
(A19) 

(A20) 
(A21) 

(A22) 



APPENDIX B: 
EVALUATION OF 

The (/3-component of the SF /' is given by 
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4 = E '^^^'^'^ 



-—I 



im 



m— — l 



All the quantities are evaluated at the position of the charge, which are assumed to be located at (r, 9^ (p). Here 



lira 



dz 



*2 



(Bl) 
(B2) 

(B3) 



and the factor AP ~ — Q"^ comes from the fact that 



dz dr 
It is easy to show that = (<i)'™)*, so is reaL Hence we can write 



= <?5(r)VM2-a2-g2 ^ Re | 

m=-l ^ 
I 

= -qS{r)y/AP-a^-Q^ E I™ 



m—-—l 



(B4) 



In Sect. U, we find that (j)^i^{z) cx P; ^(z) and cf)^2^[z) oc (z). However, the function ^{z) is related to Qi{z) 
by [H 



r(z + /i + 1) 



(z). 



Hence ^ and Qf are linearly dependent, and so we can as well choose (/)2™ oc Q; ^. In this Appendix, we set 



where {21 + 1)!! = 1 • 3 • 5 • • • (2/ + 1) and 



We have used Eq. (Elh and the fact that 



m/2 



(2; + !)!! z'+i 

3 1 



F^ 1 



r I 



(2^ + 1)!! 

21+1 



(B5) 

(B6) 
(B7) 

(B8) 
(B9) 



to obtain the second equality in Eq. (B7). Since /i ~ irwy is purely imaginary, we have — [K^)* . It follows from 
Eqs. @ and (|B7|) that 



K-^{z)^[K^{z)Y=K^{z), 



(BIO) 



Hence Ki is a real- valued function, although it is not obvious to see this from Eq. (BS). Qi{z) is then just a complex 
factor independent of z times a real- valued function, so we can choose <l)^^{z) to be real. In fact, the (fipiz) in Eq. ( |3^ ) 
is real, because it is equal to fCf(z)/z'+^. In this Appendix, however, we stick to the choice in Eq. (p7|). 
The Wronskian is given by [ESl 



W, 



1-Z2 



r (i±i£ + r 
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where we have used the identity 



r(;-^ + i) 
1 - z2 r(; + ^ + i) 



22^-5 / 1 

r(2.). r(.)r . + - 



Combining Eqs. (B6), (B7), (BIC) and (Bll), we obtain 



0'r(z)4"(z) 



1 l-z' 



[T{l + ^,+ l)pr{z)\ 



(Bll) 



(B12) 



(B13) 



Note that the quantity in th e first bracket is real, while that in the second is complex. To extract the imaginary part, 
we subtract from Eq. (B13) its complex conjugate and then divide by 2i: 



Ini 



1 I- 



z'+i (2; + i) 

1 1-^2 



V{l + ^i+ l)P^^{z) - T{1 - /i + l)Pt{z) 
2i 



r{l + fi+l)T{l- sinh(7rm7) ' 

1+1 



7r(2Z + 1)!! 



p-^{z) + -e^"'' sm{^i^T) gp(z) 

TT 



y smh{TTy) 



z'+i (2Z + 1)!! ' 
where we have used the fact that /i = imj and ||2 

_ Til + ^, + ^) ■ 

The product r(/ + /i + l)r(/ — /i + 1) can be evaluated by the identities T{x + 1) = xT{x) and |2£ 

TT 

r(i?;)r(-ij/) 

for real y. The result is 



Im 



0i'"(z)(/.!,™(2) 



an (i_ 1. 

^21 K ^2 



Substituting the above formula into Eq. (B4), we find 



^ ^ m— — l 



n(mv+/) 



[(2; + l)!! 



(2; + l)!! 



and substituting for the value of 7 we finally obtain 



7n— — l 



(2; + l)!! 



(B14) 



(B15) 



(B16) 



(B17) 



Yrr^i0,ip)Yi^{9,ip) , (B18) 



(B19) 



APPENDIX C: 



NUMERICAL EVALUATION OF THE SELF FORCE 

A very effective way to evaluate the associated Legendre functions Pj^^{z) and Qi{z) numerically is by their 
representations using hypergeometric functions. Specifically, we rewrite the associated Legendre functions as 



r(l + ^) \z-l 



z + \ 



-?,? + 1;1 + m; 



I- z 



(CI) 
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and 



In all the expressions we need to evaluate, we only have products of P^^{z) and (z) divided by their Wronskian 



determinant. Hence, we do not need to evaluate the constant factors in Eqs. ( |Cl| ) and (C2). The derivatives of 
Pj~'^{z) and Q'^{z), which are needed both for the Wronskian determinant and for the gradient of the field in the SF 
computation, can be computed using the relation 

^ 2Fi{a,b;c;z) = —2Fi{a+l,b+l;c+l;z). (C3) 



dz 



c 



It is convenient to evaluate the hypergeometric functions in Eqs. (C1-C2) using their Gauss series representation, i.e., 

2Fi(a,6;c;z) = f;M^z" |z|<l, (C4) 
n.[c)n 

where (a)n is Pochhammer's symbol. This is a useful numerical approach because of the following. For the evaluation 
of P^^(z) we need to compute a hypergeometric function for which the Gauss series is reduced to a polynomial of 
degree Hn (1 — 2;)/2, thanks to the first variable of the hypergeometric function being a negative integer. Hence, we 
only need to sum over a finite nu mber of terms, and there is no truncation error involved. We are also unrestricted by 
the radius of convergence of (p4|). This fortunate circumstance does not change for the calculation of the derivative. 
For the computation of Q^(z) we benefit from the argument being positive and smaller than unity. Consequently, in 
the Gauss series we have a sum over terms which are a (combinatorial factor) x 2;^^", where z > 1 and n is a positive 
integer. This guarantees that the series converges fast, such that in practice we do not need to sum over too many 
terms to ensure a given accuracy. In practice we evaluate all hypergeometric function to accuracy of 1 part in 10^^. 
Our computation does not become more complicated because of being imaginary. 

An alternative approach for computing Q^{z) is to use its integral representation, given by 

Qi^^)-^ 2'+ira + i)^^ li (z-t)w ^^^^ 

which is efficient numerically, because the integration interval is compact, and the integrand has pathologies neither 
inside the interval nor at its boundaries. 

The spherical harmonics can be computed either from the associated Legendre functions computed using the 
hypergeometric representation above (which is not very convenient because of the large arguments of the gamma 
functions, which are needed to be computed now), or, alternatively, from the recurrence relations for the associated 
Legendre functions as given in |]3^ . This latter approach is stable and accurate also for the very large values of Z,m 
we need (« 1000). 
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